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ABSTRACT 


L 

A  full-scale  model  of  a  double-wedge,  simulating  a  jet 
vane  used  in  guided  missile  thrust  vector  control  (TVC) ,  was 
constructed  for_ subsonic  wind  tunnel  testing  to  evaluate  the 
feasibility  of/^rtiSWin^  a  thermal  imaging  (infrared)  system 
vice  thermocouples  to  determine  a  surface  heat  transfer  coeffi¬ 
cient.  The  model  was  preheated,  then  allowed  to  cool  via 
forced  convection  to  ambient  conditions.  Surface  temperature 
readings  were  taken  using  thermocouples  and  a  thermal  imaging 
system  during  the  cooldown  process.  These  readings  were  then 
reduced  to  solve  for  a  heat  transfer  coefficient  (h) ,  using 
Newton's  Law  of  Cooling.  The  results  were  compared  with  the 
theoretical  solution  for  a  flat  plate  and  with  the  results 
obtained  by  testing  an  actual  TVC  vane  configuration  under 
the  identical  flow  conditions. 
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I .  INTRODUCTION 


A.  BACKGROUND 

The  use  of  thrust  vector  control  (TVC)  in  future  cruise 
missile  guidance  and  control  systems  is  presently  being  inves¬ 
tigated.  Removing  external  aerodynamic  steering  surfaces 
(fins)  would  allow  installation  of  a  single  weapon  configuration 
in  multiple  platforms  with  little  or  no  restructuring  of  each 
platform's  unique  launch  capabilities,  thereby  reducing  design, 
development,  and  construction  costs.  In  addition,  the  maneu¬ 
vering  envelopes  required  of  present  and  future  missiles  may 
contain  regions  in  which  external  aerodynamic  surfaces  are 
inadequate  for  control  purposes. 

The  incorporation  of  TVC  systems  for  missile  steering  re¬ 
quires  the  design  of  jet  vanes,  or  other  deflection  devices, 
that  are  mounted  in  the  rocket  exhaust  flow.  Such  devices  can 
be  rotated  to  various  angles  of  attack  in  order  to  redirect  the 
flow,  thus  maneuvering  the  missile  in  the  desired  direction. 

The  vanes  must  be  made  to  withstand  severe  steady  and  transi¬ 
ent  thermal  loads,  shock  wave  impingement,  and  erosion  from 
motor  exhaust  particulates  during  extended  use. 

One  of  the  first  steps  in  the  design  of  such  a  vane  is  to 
develop  a  model  of  the  thermal  environment  to  which  the  vane 
will  be  subjected.  This  includes  the  determination  of  surface 
heat  transfer  coef f icients ■ for  particular  vane  configurations 
under  forced  convection  conditions. 


The  following  discussion  concerns  the  development  of 
experimental  methods  to  determine  a  heat  transfer  coefficient 
on  a  heated  surface  in  subsonic  flow.  Utilizing  state-of-the- 
art  infrared  technology,  in  place  of  thermocouples,  to  determine 
surface  temperature  gradients  as  a  body  is  heated  (or  cooled, 
as  in  this  case)  would  reduce  the  time  and  effort  required  to 
develop  instriamented  models.  The  method  also  has  the  potential 
to  reduce  or  eliminate  the  need  for  elaborate  data  acquisition 
systems,  while  maintaining  an  accuracy  that  is  at  least  as 
good  as  that  obtained  from  conventional  methods.  Finally, 
the  methods  described  here  provide  a  global  as  well  as  local 
view  of  surface  heat  transfer.  Thus  it  is  possible  to  detect 
and  analyze  regions  of  particular  interest  that  might  otherwise 
go  unnoticed  in  the  data  produced  by  point  temperature 
measurements . 


B.  OBJECTIVES 

The  objectives  of  this  study  are  to: 

1)  Design  a  typical  TVC  model  instrumented  for  surface 
temperature  measurement  using  both  thermocouples 
and  thermal  imaging  technology. 

2)  Record  surface  temperatures  of  the  model  in  a  subsonic 
wind  tunnel  environment. 

3)  Compare  the  results  from  both  measurement  techniques  in 
calculating  the  heat  transfer  coefficient,  using 
Newton's  Law  of  Cooling,  and  correlating  the  results 
with  the  theoretical  solution  for  a  flat  plate  under 
identical  flow  conditions. 

4)  Utilize  the  thermal  imaging  technique  to  analyze  the 
surface  temperature  map  of  an  actual  TVC  vane 
configuration . 


5)  Present  conclusions  and  recommendations  for  future 
applications  in  this  area. 


II.  MODEL  DESIGN  AND  CONSTRUCTION 


A.  MODEL  DESIGN 

A  double-wedge  design  was  chosen  for  this  experiment  for 
several  reasons,  including  ease  of  construction,  similarity 
to  typical  TVC  configurations  used  in  previous  research,  and 
sufficient  similarity  to  a  flat  plate  to  allow  interpretation 
of  results.  A  drawing  of  the  wedge  is  included  as  Figure  1. 

The  dimensions  of  the  model  were  chosen  to  reflect  essen¬ 
tially  full-scale,  but  constrained  by  the  necessity  of  having 
a  Biot  number  less  than  0.1.  This  constraint  is  mandated  by 
the  use  of  the  lumped-heat  capacitance  analogy  in  this  experi¬ 
ment  to  determine  the  heat  transfer  coefficient  [Ref.  1]. 

This  method  of  analysis  assumes  that  the  internal  resistance 
of  the  body  is  negligible  compared  to  the  external  resistance 
and,  therefore,  a  uniform  temperature  distribution  exists 
throughout  the  model.  By  ensuring  that  the  model  has  a  Biot 
number  less  than  0.1,  this  analytical  method  will  produce 
reliable  results. 

The  convection  heat  loss  from  the  heated  model  is  assumed 
to  equal  the  decrease  in  its  internal  energy.  Thus,  the 
energy  balance  expressing  the  rate  of  heat  transfer,  q,  is: 

q  =  hA(T  -  TJ  =  ~  ^  i? 


where : 


h 


heat  transfer  coefficient,  BTUA^r-ft  -°F 


2 

A  =  surface  area,  ft 
T  i  model  temperature  at  time  t,  °F 
=  free  stream  temperature,  “F 
c  E  heat  capacity  of  model  material,  BTU/lbm-°F 
V  E  model  volume,  ft^ 

k  E  thermal  conductivity,  BTU/hr-ft-°F 

3 

p  =  density,  Ibm/ft  . 

Given  the  initial  condition  of  T  =  at  time  t  =  0,  the 
solution  to  the  differential  equation  becomes: 


T  -T^ 

^ ■  (Bl  XFO) 

where : 

Bi  Biot  number  = 

Fo  Fourier  number 

B.  CONSTRUCTION 

Copper  was  chosen  as  the  model  material  due  to  its  excel¬ 
lent  thermal  conductivity  (k  =  230  BTU/hr-f t-°F)  and  reasonably 
good  machinabil ity .  With  the  dimensions  as  shown  in  Figure  1, 

giving  V/A  =  0.01297,  the  surface  heat  transfer  coefficient 

2 

could  theoretically  rise  to  well  over  h  =  1400  ETU/hr-ft  -°F 


hV 

kA 


A,5)2t 

pc'v' 


before  the  restriction  on  Biot  number  would  be  violated.  As 
the  expected  heat  transfer  coefficient  was  less  than  60,  these 
design  decisions  ensured  excellent  accuracy  in  applying  the 
method . 

The  copper  wedge  was  instrumented  with  copper-constantin 
(type  T)  thermocouples,  30  gauge.  The  thermocouples  were 
placed  as  shown  in  Figure  1,  with  one  of  the  thermocouples 
embedded  at  depth  of  0.30  inches  to  detect  any  distinguishable 
temperature  difference  between  the  surface  and  the  interior. 

As  can -be  seen  in  the  subsequent  tables  and  graphs,  the  assump¬ 
tion  of  Newtonian  cooling  proved  to  be  accurate.  The  wedge 
was  also  covered  with  a  thin  coat  of  flat  black  paint  to  in¬ 
crease  the  emissivity  of  the  surface  to  a  minimum  of  0.95 
[Ref.  1]  in  order  to  enhance  the  thermal  image  detection  of  the 
infrared  system. 

The  thermocouples  were  attached  to  a  dual-channel  strip 
chart  recorder,  and  to  a  digital  temperature  sensor.  The 
digital  readout  was  used  to  obtain  reference  values  of  wedge 
temperature  during  heatup,  and  free  stream  temperatures  during 
testing.  All  instruments  and  thermocouples  were  calibrated 
from  32°F  to  150°F  using  the  Department  of  Mechanical  Engineer¬ 
ing  calibration  facility. 

The  wedge  was  then  mounted  on  two  aluminum  stanchions, 
attached  to  the  plexiglass  insert  in  the  floor  of  the  wind 
tunnel  test  section.  Within  the  plexiglass,  a  six-inch  diam¬ 
eter  hole  was  cut  for  insertion  of  the  infrared-transparent 


window  through  which  the  thermal  imaging  system  could  view 
the  wedge  surface.  This  arrangement  is  shown  in  Figure  2. 
Plexiglass  sidewalls  were  added  to  form  a  thermal  barrier  to 
attempt  to  negate  any  spanwise  heat  transfer  from  the  wedge 
sides,  and  to  form  a  channel  for  uniform  air  flow  past  the 
wedge.  The  model  was  heated  by  using  an  electric  pad,  capable 
of  raising  temperatures  from  ambient  to  approximately  140°F 
within  a  thirty  minute  time  frame. 


III.  THERMAL  IMAGING  SYSTEM 


A.  DESCRIPTION 

The  thermal  imaging  system  used  in  this  procedure  is  a 
Probeye  Thermal  Video  System,  Series  4300,  manufactured  by 
Hughes  Aircraft  Company  [Ref-  2] .  This  system  provides  real¬ 
time  temperature  maps  of  objects  through  16  color  bands  dis¬ 
played  on  a  color  monitor.  The  temperature  range  of  the  system 
is  from  -4°F  to  536‘’F  (-40°C  to  +280°C)  .  Each  color  on  the 
scale  can  represent  temperature  sensitivities  from  0.5  to 
36°F  (0.5  to  20°C) ,  depending  upon  the  desired  overall  tempera¬ 
ture  range.  Emissivity  adjustments  can  be  made  from  0.0  to 
1.0  in  increments  of  0.01,  in  the  spectral  wavelength  band  of 
2.0  to  5.6  umeters.  The  unit  includes  a  movable  cursor  which 
can  be  placed  at  any  location  on  the  viewed  surface,  giving  a 
digital  readout  of  the  temperature  at  that  location,  and  a 
real  time  clock/calendar  precise  to  0.1  seconds  for  determining' 
time-variant  conditions. 

The  use  of  this  imaging  system  also  required  the  use  of 
an  infrared  optic  window  for  viewing  the  model  while  maintain¬ 
ing  tunnel  integrity.  A  six-inch  diameter,  two  millimeter 
thick  window  fabricated  from  magnesium  fluoride  was  obtained 
from  Eastman  Kodak  Company.  This  type  of  material,  tradename 
IRTRAN  1,  was  chosen  to  match  the  optical  wavelength  charac¬ 
teristics  of  the  imaging  system.  As  manufactured,  the  window 


allowed  over  95%  transmittance  in  the  2.0  to  5.8  ymeter 
range,  while  capable  of  withstanding  pressure  differentials 
of  1  atmosphere  at  25°C,  with  a  factor  of  safety  of  four 
[Ref .  3]  . 
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IV.  EXPERIMENTAL  PROCEDURE 


A.  WIND  TUNNEL  TESTING 

The  model  assembly,  as  shown  in  Figure  2,  was  installed 
in  the  test  section  of  the  subsonic  wind  tunnel.  This  wind 
tunnel  is  a  pull-down  system  capable  of  producing  test  section 
free  stream  velocities  of  up  to  335  feet/second  by  means  of  a 
variable  speed  fan  located  at  the  tunnel  exit.  The  free-stream 
turbulence  level  of  the  tunnel  (0.2%-0.5%)  was  deemed  suita¬ 
ble  for  the  purposes  of  this  experiment  [Ref.  4]. 

At  the  start  of  each  test,  with  the  tunnel  fan  off,  the 
wedge  was  wrapped  with  the  heating  pad,  and  brought  up  to 
temperatures  of  approximately  130“F.  The  pad  was  then  removed, 
wind  tunnel  access  doors  closed,  and  the  fan  started.  Runs 
were  made  at  50%  and  75%  fan  speed,  corresponding  to  200  and 
306  feet/second,  respectively.  Free  stream  velocities  were 
calculated  using  the  installed  manometer,  calibrated  in  inches 
H2O.  Free  stream  temperatures  were  determined  by  both  a 
calibrated  thermometer  and  a  thermocouple  inserted  in  the  flow. 

The  cooling  process  was  then  recorded  using  the  installed 
thermocouple  output  to  the  strip  chart  recorders,  and  the  ther¬ 
mal  imaging  system  connected  to  a  video  tape  recorder  for 
simultaneous  data  acquisition.  in  initial  tests,  the  tempera¬ 
ture  sensitivity  settings  on  the  imaging  system  were  selected 
at  5.0°F,  which  yielded  an  80  degree  range  from  50°F  to  llO^F, 


to  view  the  entire  cooling  process.  In  subsequent  tests,  a 
sensitivity  of  2.0“F  was  selected,  giving  a  32  degree  range 
from  90°F  to  122°F.  This  lower  sensitivity/range  setting  was 
necessary  to  obtain  a  good  agreement  with  the  thermocouple 
data . 

After  a  considerable  amount  of  system  testing,  this  process 
was  undertaken  four  times  at  each  of  the  two  fan  speeds  for 
data  evaluation  and  correlation. 

B.  DATA  REDUCTION 

Applying  the  formula  for  the  lumped-capacitance  analysis, 
the  transient  data  from  the  thermocouple  readings  and  the 
thermal  imaging  system  were  used  to  deduce  the  corresponding 
heat  transfer  coefficients.  The  calculations  were  performed 
using  TK]  SOLVER,  a  microcomputer  spreadsheet  [Ref.  5].  Tables 
1  and  2  exhibit  the  program  equations  and  variables,  which 
provide  not  only  the  heat  transfer  coefficient,  but  also  an 
estimate  of  the  uncertainty  of  that  calculation,  and  the 
corresponding  Biot  number,  for  each  of  the  two  sensor  types. 

The  uncertainty  analysis  of  the  heat  transfer  coefficient 
calculation  is  included  as  Appendix  A. 

The  reduction  of  the  thermocouple  data  was  fairly  straight¬ 
forward,  as  the  output  was  temperature  versus  time.  However, 
reducing  the  thermal  imaging  system  data  required  running  the 
video  tape  recording  frame  by  frame  (each  frame  corresponding 
to  1/60  second)  to  determine  the  temperature  at  a  given  time. 


This  method  generated  the  largest  uncertainty  in  the  calculation 
because  the  visual  color  bands  and  digital  readout  of  the  sys¬ 
tem  changed  only  when  the  temperature  of  the  body  changed 
more  than  the  sensitivity  setting.  For  example,  at  a  5.0°F 
setting,  the  readout  would  persist  for  perhaps  three  seconds 
of  test  time  before  changing  from  110°F  to  105°F,  while  in 
reality,  as  displayed  by  the  thermocouple  readings,  the  surface 
temperature  could  be  any  temperature  within  this  range  at  a' 
particular  instant.  (This  problem  is  analogous  to  that  due 
to  a  large  interval  size  in  a  finite-difference  calculation.) 
There  also  existed  occasions  when,  particularly  at  temperatures 
near  thermal  equilibrium,  the  readout  would  take  relatively 
long  periods  of  time  to  stabilize,  which  made  determining  the 
exact  temperature/time  data  difficult.  However,  these  diffi¬ 
culties  were  overcome  by  operator  experience  and  use  of  the 
narrower  sensitivity  setting  for  correlating  data  in  smaller 
segments  rather  than  over  the  entire  cooling  span. 


I 

I  V.  EXPERIMENTAL  RESULTS 

A.  FLAT  PLATE  CORRELATION 

The  analysis  of  a  flat  plate  of  identical  material,  length, 
and  width  in  identical  flow  conditions  was  used  to  attain  a 
degree  of  confidence  in  the  accuracy  of  the  results  of  the 
wedge  experiment.  This  analysis  is  based  upon  determining 
^  the  Reynolds  and  Prandtl  numbers  for  the  flow,  and  solving  the 

following  equation  for  the  average  Nusselt  number:  [Ref.  6] 

I  Nu  =  0.664  Re^^^Pr^/^ 

The  Nusselt  number  was  then  used  to  determine  the  average 
^  heat  transfer  coefficient  from: 

r 

Nu  k 

h  =  Jl 

where 


I  =  length  of  flat  plate,  ft. 

These  procedures  resulted  in  a  heat  transfer  coefficient  of 

22  BTU/hr-f t^- ‘’F  for  free  stream  velocity  of  200  ft/sec,  and 
2 

28  BTU/hr-ft  -°F  at  velocity  of  306  ft/sec.  The  results  of 


this  calculation  are  shown  in  Appendix  B. 


B .  THERMOCOUPLE  RESULTS 

Using  the  strip  chart  recorder  traces,  the  heat  transfer 

coefficient  at  the  surface  thermocouple  location  was  determined 

from  the  lumped-capacitance  solution  of  exponential  temperature 

decrease  with  time.  For  the  free  stream  velocity  of  200  ft/ 

2 

sec,  h  varied  from  27.5  to  33.2  BTU/hr-ft  -°F,  with  an  average 

2 

uncertainty  band  of  ±1.75  BTU/hr-ft  -°F.  At  velocity  of 

306  ft/sec,  h  had  values  ranging  from  41.03  to  42.7  BTU/ 

2  2 
hr-ft  -®F,  with  average  uncertainty  of  ±1.80  BTU/hr-ft  -°F. 

These  calculations  are  shown  in  Tables  4-11. 

C.  THERMAL  IMAGING  SYSTEM  RESULTS 

As  shown  in  Tables  12-19,  using  the  same  analysis  and 

temperature  range,  h  values  were  determined  to  be  27.7  to 

33.5  BTU/hr-ft^-°F,  with  an  uncertainty  band  of  ±1.9  BTU/ 

2 

hr-ft  -°F  for  the  lower  free  stream  velocity,  while  values  of 
39.9  to  41.8  were  obtained  at  higher  speed.  Figures  3  and  4 
show  the  comparison  of  wedge  surface  temperature  versus  time 
for  the  two  sensors  at  each  velocity,  and  Figures  5  and  6 
show  the  calculated  values  from  each  sensor  versus  time,  again 
at  each  speed. 

It  is  important  to  note  that  this  measurement  corresponds 
to  conditions  at  a  particular  point  on  the  wedge,  that  is,  the 
location  of  the  corresponding  thermocouple.  Further  testing, 
with  the  cursor  positioned  at  other  points,  would  allow  the 
development  of  a  complete  heat  transfer  map  of  the  wedge  sur¬ 
face.  However,  viewing  of  the  entire  wedge  surface  showed  that 


the  transient  isotherm  variation  was  essentially  uniform  over 
a  large  portion  of  the  surface  (with  the  exception  of  narrow 
bands  in  the  three-dimensional  boundary  layer  regions  near  the 
sidewalls) .  Figure  7  is  a  still  photograph  of  the  wedge  taken 
from  the  video  output,  which  clearly  exhibits  this  condition. 
Thus  it  may  be  deduced  that  the  values  determined  above,  at 
the  thermocouple  location,  are  representative  of  conditions 
over  a  major  portion  of  the  wedge  surface.  The  ability  to  make 
such  observations  is  a  significant  virtue  of  the  infrared 
thermography  method. 

D.  TVC  VANE  ANALYSIS 

An  actual  TVC  vane  was  received  from  NWC  China  Lake,  Ca 
for  analyzing  the  effect  of  complex  geometry  on  convective  heat 
transfer  using  the  thermal  imaging  system. 

This  vane,  shown  in  Figure  8,  was  constructed  of  copper- 
impregnated  tungsten  (80%  W,  20%  Cu) ,  with  the  following 
material  properties: 


p  =  1076.544 


148  _  c 

®  hr-ft-°F 


0.037 


BTU 

lbm-°F 


The  vane  volume  was  calculated  to  be  2.295  xlO  ^  ft^,  with 

-1  2 

a  corresponding  surface  area  of  1.5551  x lo  ft  .  All  values 
were  then  entered  into  a  formula  sheet  as  before,  shown  in 
Table  3. 


The  vane  was  mounted  on  two  aluminum  stanchions  with  the 
same  baseplate  used  by  China  Lake  in  their  experiments.  The 


thermal  imaging  system  was  set  up  to  view  the  vane  through  the 
side  port  of  the  subsonic  wind  tunnel  with  the  IRTRAN  optic 
window  installed.  The  identical  heating  procedure  was  used, 
however,  for  reasons  unknown,  the  pad  would  only  reach  105 °F. 
This  was  deemed  sufficient,  as  it  still  gave  a  40  degree  tem¬ 
perature  differential.  The  wind  tunnel  was  operated  at  50% 
flow  speed  (200  ft/sec). 

Three  test  runs  were  performed,  with  the  cursor  at  the 
midspan  for  runs  1  and  3,  and  on  the  leading  edge  at  the  base¬ 
plate  boundary  for  run  2 .  The  cursor  was  moved  for  run  2  to 
investigate  and  analyze  a  more  rapid  temperature  gradient 

noticed  in  that  area  in  the  previous  run.  Run  1  yielded  an 

2 

average  h  of  41.5  BTU/hr-ft  -°F  +3.23,  run  2  yielded  an  average 
of  54.8  +4.41,  and  run  3  gave  an  average  of  42.4  ±3.9.  The 
largest  Biot  number  was  calculated  as  0.0063,  well  within  the 
lumped-capacitance  analysis  criteria. 

Figure  9  is  a  color  thermograph  of  the  vane  during  cool¬ 
down.  As  compared  to  the  thermogram  of  the  double  wedge 
(Figure  7),  the  temperature  map  is  fairly  uniform  over  the 
body,  with  the  exception  of  faster  cooling  at  the  leading  and 
trailing  edges,  and  an  interesting  area  near  the  top  of  the 
vane  below  the  apex.  While  vane  thickness  would  explain  the 
shape  of  most  of  the  isotherms,  the  upper  region  undoubtedly 
is  affected  by  a  unique  flow  field,  which  increases  the  heat 
transfer  coefficient  by  approximately  8%  (Table  22).  Tables 
20-22  show  the  results  of  these  tests,  and  Figures  10  and  11 
are  graphs  of  temperature  and  h  versus  time,  respectively. 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  results  obtained  from  the  thermal  imaging  system  are, 
within  a  small  experimental  uncertainty,  in  agreement  with 
the  results  of  the  thermocouples  for  all  conditions  tested 
here.  Both  experimental  method's  results  were  also  suffi¬ 
ciently  close  to  the  theoretical  flat  plate  solution  to  allow 
confidence  in  their  accuracy,  primarily  due  to  the  small  wedge 
angle  used. 

The  primary  disadvantage  of  the  thermal  imaging  system  is 
the  increase  in  size  of  the  uncertainty  band  as  the  tempera¬ 
ture  sensitivity  is  widened.  For  large  sensitivities  of  5 
degrees  or  more,  the  results  become  comparatively  less  accurate, 
again  similar  to  using  large  step  sizes  in  finite  element 
calculations.  Future  users  of  this  system  in  experiments 
similar  to  this  one  will  have  to  adjust  the  sensitivity  settings 
with  respect  to  the  overall  temperature  range  desired,  most 
likely  by  trial  and  error,  to  get  the  degree  of  accuracy 
required . 

Nevertheless,  the  results  obtained  in  this  experiment 
lead  to  the  conclusion  that  the  thermal  imaging  system,  when 
properly  employed,  would  be  a  very  good  alternative  to  thermo¬ 
couples  in  applications  similar  to  this  one.  It  provides 
accurate  local  and  global  readings  with  no  need  for 


complicated  data  acquisition/reduction  equipment  and  extensive 
instrumentation.  The  system  is  transportable,  very  reliable, 
and  easy  to  operate.  Although  the  system  requires  special 
infrared  optic  material  for  viewing  models  in  an  enclosed 
environment,  several  vendors  were  found  who  are  capable  of 
providing  the  material  in  various  shapes  and  sizes  at  a  reason 
able  expense  in  time  and  money. 

B.  RECOMMENDATIONS 

The  findings  of  this  research  lead  to  a  recommendation  for 
using  this  system  for  surface  temperature  measurements  in  the 
next  phase  of  this  project,  using  the  wedge  model  in  a  super¬ 
sonic  wind  tunnel  for  the  same  heat  transfer  coefficient  calcu 
lation.  This  model  was  designed  for  ease  of  installation  in 
the  supersonic  tunnel,  when  it  becomes  available,  and  an  optic 
window  has  been  obtained  that  will  allow  infrared  viewing  of 
the  model  while  withstanding  the  environment  (temperature  and 
pressure  differentials)  in  supersonic  flow.  Particular  problems 
that  might  be  encountered  in  supersonic  applications  include: 

1)  Sufficient  difference  in  model  and  free  stream 
temperatures  to  cause  adequate  heat  transfer — the 
model  may  have  to  be  heated  prior  to  testing  by 
some  internal  means  (installation  of  wires  for 
Joulean  heating)  to  increase  the  recovery  factor. 

2)  The  on-site  supersonic  tunnel  has  a  significantly 
smaller  run-time  available  due  to  its  pressure 
vessel  configuration,  thus  yielding  short  run 
durations  (a  few  minutes)  and  long  (hours  vice 
minutes)  pressurization  durations  between  tests. 

3)  As  other  model  shapes  and  materials  are  considered, 
the  necessity  for  maintaining  the  Biot  number  less 
than  0.1  could  become  critical,  or  another  method 


of  analysis  would  have  to  be  chosen  if  the  lumped- 
capacitance  assumption  is  not  valid. 

It  is  also  recommended  that  the  next  experimenter  familiar 
ize  himself  with  the  system  thoroughly,  and  invest  time  in 
viewing  the  video  tape  outputs  of  the  system  to  get  experience 
in  temperature  determination.  This  experience  will  prove  to 
be  a  great  factor  in  enhancing  the  accuracy  of  the  results. 


APPENDIX  A 


UNCERTAINTY  ANALYSIS  CALCULATIONS 


T  -T 


0  = 


T  _T 

o  °° 


,hAt . 


Formula  used  in  calculation 


h  = 


^ A^t^  ^  solving  for  h 


h  =  F(0,t,A,V,p,c) 


Ah  = 


r,9h  ,-.2  ^  ,3h  ..,2  ^  .9h  ...2  ,9h  .,,.2 


1/2 


where  A's  are  the  uncertainty  of  the  measured  quantities. 


(1)  gg  A0 


££|[(-  ^  +  (b  an)  2, 


1/2 


where ; 


N  = 


T  -T  ,  D  =  T  -T  ,  AN  =  (AT^  +  AT^) 

00  '  O  oo'  CO 


1/2 


AD  =  (AT^  +  AT^ 
o  °0 


1/2 


,  AT  =  1.0  °F,  AT  =  1.0  °F 


AT  =  2.0  °F  for  IR  data,  AT  =  1 . 0  ®F  for  thermocouple 


2.7778  E  -5  hours 


(2) 


5t  “ 


In  0  pcV 

—  "  ^  a  u  ^ 

At 


where : 

At  =  0.1  seconds  = 

from  IR  system  timer  precision. 

(3)  !§■  AA  =  -  Q  P-gJi  aa 

3A 

where : 

AA  =  ■  0.001  in^  =  6.9444E-6ft^ 

from  machining  precision. 

(4)  111  AV  =  -  P,.2 

(‘i)  3v  At  ' 


where : 


AV  =  l.OE-9  in^  =  5.7878E-13  ft^ 


from  machining  precision. 


,  r-  .  9h  A 

(5)  ^  Ap  = 


«,n  0c  V 
At 


Ap  , 


where : 


Ap 


6.9082 


Ibm 


[Ref.  7] 


(6) 


3h 


Ac 


Jin  0  pV 

A  t 


AC  , 


where : 


Ac 


0.001 


BTU 

lbm-°F 


[Ref . 


7] 


All  above  equations  were  included  in  the  data  reduction 
programs  to  solve  for  the  uncertainty  in  h  (delh) . 
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APPENDIX  B 


FLAT  PLATE  SOLUTION 

Given:  T  =  60°F  =  520°R  1"H„0  =  0.0361  psi 

OO  2  ^ 

Find:  h  from  the  empirical  flat  plate  in  parallel  flow 

equation: 


Nu  =  0.664  Pr^^^ 

X 


h  = 


Nu  k 


[Ref.  6] 


(1)  Re  :  For  T  =  Sao^R:  [Ref.  7] 

X  <» 


y  .  =  1.2039E-5  Pr 

^air  ft-sec 


=  0.71  k  =  0.01561 


hr- 


Y  .  =  0.0743 

'  air 


Ibm 


ft' 


X  =  2.50  in  =  0.2083  ft 


y2  ^  AP_25  ^ 


(8.9"H20)  (0.0361)  (64.348)  (144] 
070743 


U  =  40080.4 


ft' 


2  ' 


=  200.2  —  for  experimental 

low  speed  flow 


f  t 

306.1  — 
s 


BTU 

ft-°F 


For  Ap 


20.8"  H-,0,  U 

2  OO 


for  experimental 
high  speed  flow 


yU  X 

•  m 


(0.0743)  (200.2)  (0.2083) 
1.2039E-5 


=  2.573E5 


(0.0743) (306.1) (0.2083) 
1.2039E-5 


=  3.935E5 


Assuming  =  5.0E5,  laminar  flow  conditions  exist  in  both 

cases . 


(2)  Nu:  For  low  speed  flow: 


Nu  =  0. 664  (2. 574E5)^'^^(0. 71)^/^)  =  300.5 


For  high  speed  flow: 


Nu  =  0. 664  (3. 935E5)-"^"-(0. 71)-"^  =  371.6 


1/3  . 


(3)  h:  For  low  speed  flow: 


Nu  k  _  (300.5)  (0.01561) 

X  0.2083 


=  22.5 


hr-ft  -°F 


For  high  speed  flow: 


h  = 


Nu  k  _  (371.6)  (0.01561 

X  0.2083 


■^  =  27.8 


hr-ft  -°F 


Vv*  ■  ■  \  •  ■- 
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APPENDIX  C 


TABLES 


TABLE  1 

FORMULA  SHEET  FOR  THERMOCOUPLE  DATA  REDUCTION 


St 

Input 

Name 

Output  Unit 

Comment 

L 

theta 

dimensionless 

temperature 

L 

Tt 

F 

temperature  at  time  t 

Tinf 

F 

free  stream 
temperature 

Ti 

F 

initial  wedge 
temperature 

T 

h 

btu 

heat  transfer 
coefficient 

JLi 

hr-ft^-°F 

L 

delh 

a 

btu 

uncertainty  in  h  value 
wedge  wetted  perimeter 

r 

.27777778 

hr-ft2-<>F 

ft"2 

L 

t 

hours 

558 

ro 

lbm/ft^3 

copper  density 

.091 

c 

btu/lb-F 

specific  heat  capacity 

.00346013 

V 

ft^3 

wedge  volume 

L 

q 

btu/hr 

heat  transfer 

230 

k 

btu 

hr-ft-°F 

wedge  conductivity 

L 

BIOT 

L 

time 

seconds 

u 

(dh/dtheta) * (delta  theta) 

vv 

(dh/dt) * (delta  t) 

w 

(dh/dA) * (delta  a) 

X 

(dh/dV) * (delta  V) 

y 

(dh/dro) * (delta  ro) 

z 

(dh/dc )* (delta  c) 

TABLE  1  (CONTINUED) 


S  Rule 


*  theta 


(Tt-Tinf ) 
(Ti-Tinf ) 


*  theta 


exp 


(-  (h*a*t) 
(ro*c*V) ) 


*  q  = 

*  BIOT 

*  time 


h*a* (Tt-Tinf) 

h*V 

(k*a) 

=  t*3600 


*  u  =  (ro*c*V)/(theta*at))*((1.4142*(Tt-Tinf)/(Ti-Tinf)"2) 

+  (1.8868/ (Ti-Tinf )) ''2) 

*vv  =  (2.777777e-5*ln(theta)*ro*c*V/(a*t"'2)  ) 

*  w  =  (In  (theta)  *ro*c*V*6 .9444e-6/ (a''2*t)  ) 

*  X  =  (In (theta) *ro*c*5. 787e-13/ (a*t) ) 

*  y  =  (In (theta) *c*V*6 .9082/ (a*t) ) 

*  z  =  (In(theta) *V*ro*.001/(a*t) ) 

*  delh  =  sort  (u''2+vv^2+w''2+x''2+y"2  +  z^2) 


TABLE  2 


FORMULA  SHEET  FOR  INFRARED  SYSTEM  DATA  REDUCTION 


St  Input 


27777778 


.091 

.00346013 


S  Rule 


Name 

theta 


Tinf 


Output  Unit 


delh 


BIOT 

time 


theta  = 


(Tt-Tinf ) 
(Ti-Tinf ) 


*  uu  h*a*t 

*  theta  =  exp(T - ^ 

^  (ro*c*V) 


*  q  =  h*a* (Tt-Tinf ) 


hr-ft  -“F 
btu 

hr-ft^-°F 

ft''2 

hours 

lbm/ft^3 

btu/lb-F 

ft^3 

btu/hr- 

btu 

hr-ft-°F 

seconds 


Comment 

dimensionless  temperature 

temperature  at  time  t 

free  stream  temperature 

initial  wedge  temperature 

heat  transfer  coeffi¬ 
cient 

uncertainty  in  h  value 
wedge  wetted  perimeter 

copper  density 
specific  heat  capacity 
wedge  volume 
heat  transfer 
wedge  conductivity 


(dh/dtheta) * (delta  theta) 
(dh/dt) * (delta  t) 

(dh/dA) * (delta  a) 

(dh/dV) * (delta  V) 

(dh/dro) * (delta  ro) 
(dh/dc) * (delta  c) 


TABLE  2  (CONTINUED) 


BIOT 

time 

u  = 

vv  = 

w  = 

X  = 

y  = 

2  * 


h*V 

(k*a) 

=  t*3600 

(ro*c*V/(theta*a*t)  )*{  (1.4142*  (Tt-Tinf)  /  (Ti-Tinf)  ''2) 
+  (2. 2361/ (Ti-Tinf) )^2) 

(2.777777e-5*ln(theta) *ro*c*V/ (a*t"2) ) 

(In  (theta )  *ro*c*V*6 .9444e-6/(a''2*t) ) 

(In (theta) *ro*c*5 . 787e-13/ (a*t) ) 

(In(theta) *c*V*6 .9082/ (a*t) ) 

(In  (theta) *V*ro* .001/ (a*t) ) 


delh  =  sqrt  (u''2+vv^2+w^2+x^2+y''2+z''2) 


TABLE  3 


FORMULA  SHEET  FOR  TVC  VANE  DATA  REDUCTION 


S  Rule 
*  theta 


(Tt-Tinf ) 
(Ti-Tinf ) 


* 


theta  =  exp ( 


-(h*a*t) 

(ro*c*V) 


) 


*  q  =  h*a* (Tt-Tinf ) 


*  BIOT 


h*V 


*  time  =  t*3600 


*  u  =  {ro*c*V/ (theta*a*t) )*( (1.4142* (Tt-Tinf )/ (Ti-Tinf ) "2) 
+  (2.2361/(Ti-Tinf) )"2) 


*  vv  =  (2.777777e-5*ln(theta)*ro*c*V/(a*t^2) ) 

*  w  =  (In  (theta)  *ro*c*V*6.9444e-6/ (a''2*t)  ) 

*  X  =  (In (theta) *ro*c*5 .787e-13/ (a*t) ) 

*  y  =  (In (theta) *c*V*6 . 9082/ (a*t) ) 

*  z  =  (In(theta) *V*ro* .001/ (a*t) ) 


*  delh 


sqrt  (u"2+vv^2+w''2+x''2+y''2+z"2) 


TABLE  3  (CONTINUED) 


St 

Input 

Name 

Output  Unit 

L 

theta 

L 

81 

Tt 

F 

61 

Tinf 

F 

97 

Ti 

F 

L 

h 

btu 

L 

delh 

hr-ft  -»F 

.155507 

a 

ft^2 

L 

t 

hours 

1076.544 

ro 

lbm/ft^3 

.037 

c 

btu/lb-F 

. 00229545 

V 

ft^3 

L 

L 

148 

q 

k 

BIOT 

btu/hr 

btu 

hr-ft-°F 

time 


Comment 

dimensionless  temperature 

temperature  at  time  t 

free  stream  temperature 

initial  TVC  vane 
temperature 

heat  transfer  coefficient 
uncertainty  in  h  value 
TVC  vane  surface  area 

TVC  vane  density-80%  W, 
20%  Cu 

specific  heat  capacity 
TVC  vane  volume 
heat  transferred 
TVC  vane  conductivity 


L 

L  47.3 


seconds 


TABLE  4 


RESULTS  OF  THERMOCOUPLE  CALCULATIONS,  LOW  TUNNEL  SPEED 


Run : 


#1 


Tunnel  Air  Velocity: 


201.3 


ft 


sec 


Average  Biot  Number: 


0.00164 


Average  Uncertainty  in  h:  ±2.28 


btu 


hr-ft^-°F 


Time  (Seconds) 


Temperature  (°F) 


h  ( 


btu 


hr-ft-“F 


0.0 

117.5 

6.7 

112.5 

14.5 

107.5 

22.7 

102.5 

32.0 

97.5 

43.0 

92.5 

30.92 

30.00 

30.32 

30.42 

30.21 


55.0 


87.5 


30.54 


TABLE  5 


RESULTS  OF  THERMOCOUPLE  CALCULATIONS,  LOW  TUNNEL  SPEED 


Run : 


#2 


Tunnel  Air  Velocity: 


193.3 


ft 

sec 


Average  Biot  Number: 


0.00156 


Average  Uncertainty  in  h:  ±1.95 


btu 


hr-ft  -°F 


V  u  /  btu  . 

Time  (Seconds)  Temperature  (°F)  n  t-op 


0.0 

114.5 

3.5 

112.5 

24.32 

5.6 

110.5 

30.99 

9.3 

108.5 

28.56 

12.8 

106.5 

28.24 

16.0 

104.5 

28.85 

21.0 

102.5 

26.97 

24.6 

100.5 

27.48 

I 
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TABLE  6 


RESULTS  OF  THERMOCOUPLE  CALCULATION,  LOW  TUNNEL  SPEED 


Run  : 


#3 


Tunnel  Air  Velocity; 


197.6 


ft 

sec 


Average  Biot  Number; 


0.00159 


Average  Uncertainty  in  h:  ±1.53 


btu 


hr-ft  -°F 


Time  (Seconds) 

0.0 
10 . 0 
20.0 
30.0 
40.0 
50.0 


Temperature  ( ° F ) 

120.0 

113.3 

107.0 

101.7 

96.7 

92.8 


h 


,  btu 
^hr-ft-” 


28.70 
29.47 
29 . 39 
29.92 
29 . 54 


60 . 0 


89.2 


29 . 54 


TABLE  7 


RESULTS  OF  THERMOCOUPLE  CALCULATIONS,  LOW  TUNNEL  SPEED 


Run : 

Tunnel  Air  Velocity: 
Average  Biot  Number: 
Average  Uncertainty  in  h: 


#4 

207.9 

sec 

0.00182 


±  2.41 


btu 

hr-ft^-°F 


Time  (Seconds) 

0.0 

5.0 

10.0 

15.0 

20.0 

25.0 

30.0 

35.0 

37.0 


Temperature  (°F) 

111.0 

107.7 

104.3 

101.4 

98.7 
96.2 

93.7 
91.6 
91.0 


h  ( 


btu 

hr-ft-“F 


32.44 

34.23 

33.72 

33.71 

33.69 

33.82 

33.57 

33.17 


40 


TABLE  8 


RESULTS  OF  THERMOCOUPLE  CALCULATIONS,  HIGH  TUNNEL  SPEED 


Run: 

Tunnel  Air  Velocity: 
Average  Biot  Number : 
Average  Uncertainty  in  h : 


#5 

308.2  — 
sec 

0.0022 


±  1.88 


btu 

hr-ft^-°F 


Time  (Seconds)  Temperature  (°F)  h  (hr-ft-°F 


0.0  119.0 


5.0 

114.3 

40.37 

10.0 

109.7 

41.45 

: 

15.0 

105.5 

41.73 

20.0 

101.7 

41.95 

ini 

25.0 

98.3 

42.03 

30.0 

95.2 

42.00 

‘Sv 

•  I 

35.0 

92.2 

42.26 

•* 

37.0 


91.0 


42.66 


TABLE  9 

RESULTS  OF  THERMOCOUPLE  CALCULATIONS,  HIGH  TUNNEL  SPEED 


Run : 

#6 

Tunnel  Air  Velocity: 

303.8  — 
sec 

Average 

Biot  Number: 

0.0022 

Average 

Uncertainty 

in  h:  ±2.11  - * - 

hr-ft  -°F 

Time  (Seconds) 

Temperature  {°F)  h 

,  btu 
^hr-ft- 

0.0 

119.0 

5.0 

114.5 

38.15 

10.0 

110.3 

38.64 

15.0 

106.7 

37.48 

20.0 

102.1 

40.84 

25.0 

98.5 

41.52 

30.0 

95.6 

41.07 

35.0 

92.4 

41.88 

38.0 


91.0 


41.54 


TABLE  10 


RESULTS  OF  THERMOCOUPLE  CALCULATIONS,  HIGH  TUNNEL  SPEED 


Run: 


Tunnel  Air  Velocity: 


Average  Biot  Number: 


Time  (Seconds] 


#7 


306.0 


ft 


sec 


0.0022 


Average  Uncertainty  in  h:  ±2.28 


btu 


hr-ft  -°F 


Temperature  (°F) 


h  ( 


btu 


hr-ft-°F 


TABLE  11 


RESULTS  OF  THERMOCOUPLE  CALCULATIONS,  HIGH  TUNNEL 


Run : 


#8 


Tunnel  Air  Velocity: 


306.0 


ft 

sec 


Average  Biot  Number: 


0.0022 


Average  Uncertainty  in  h:  ±1.96 


btu 


hr-ft  -°F 


Temperature  (°F) 


h  { 


hr- 


SPEED 


btu 

ft-°F 


Time  (Seconds) 


TABLE  13 


RESULTS  OF  INFRARED  SYSTEM  CALCULATIONS,  LOW  TUNNEL  SPEED 


Run : 

Tunnel  Air  Velocity: 
Average  Biot  Number: 
Average  Uncertainty  in  h: 


Time  (Seconds)  Temperature  (°F)  h 


0.0 

114.5 

3.1 

112.5 

28  .40 

6.1 

110.5 

28.46 

9.5 

108 .5 

27  .96 

13.3 

106.5 

27  .18 

17 .0 

104  .  5 

27  .15 

21.5 

102  .5 

27  .33 

#2 


193.3 


ft 

sec 


0.00151 


±  2.16 


btu 


hr-ft  -°F 


24.4 


100.5 


27.71 


TABLE  14 


RESULTS  OF  INFRARED  SYSTEM  CALCULATIONS,  LOW  TUNNEL  SPEED 

Run;  #3 

Tunnel  Air  Velocity:  197.6 

^  sec 

Average  Biot  Number:  0.00161 

Average  Uncertainty  in  h:  ±1.98  - - 

hr-ft"^-°F 


jerature 


h  '  ^ 

^hr-ft-°F^ 


Time  (Seconds ) 


TABLE  15 


RESULTS  OF  INFRARED  SYSTEM  CALCULATIONS,  LOW  TUNNEL  SPEED 

Run:  #4 

Tunnel  Air  Velocity:  207.9 

Average  Biot  Number:  0.0018 

Average  Uncertainty  in  h:  ±2.14  - - 

hr-ft  -®F 


Time  (Seconds) 


Temperature  (°F) 


h  r  \ 

^hr-ft-°F^ 


109.0 

34.61 

107.0 

33.02 

105.0 

32.35 

103.0 

33.48 

101.0 

32.64 

34.44 

34.19 

33.65 


33.53 


TABLE  16 


RESULTS  OF  INFRARED  SYSTEM  CALCULATIONS,  HIGH  TUNNEL  SPEED 


Run; 

Tunnel  Air  Velocity; 
Average  Biot  Number : 
Average  Uncertainty  in  h; 


#5 


308.2 


ft 

sec 


0.0021 


±  1.85 


btu 

hr-ft^-°F 


Time  (Seconds)  Temperature  (°F) 


h 


,  btu 
^hr-ft-°F 


) 


0.0 

119.0 

4.7 

115.0 

35.91 

8.9 

111.0 

39.44 

13.7 

107.0 

40.08 

19.3 

103.0 

39.70 

23.8 

99.0 

42.27 

29.6 

95.0 

43.05 

91.0 


41.75 


TABLE  17 


Time  (Seconds) 


TABLE  18 


RESULTS  OF  INFRARED  SYSTEM  CALCULATIONS,  HIGH 

Run : 

Tunnel  Air  Velocity: 

Average  Biot  Number: 

Average  Uncertainty  in  h: 

Time  (Seconds)  Temperature  (°F) 


#7 


306.0 


ft 

sec 


0.0021 


±  1.86 


btu 


hr-ft^-®F 


0.0 

119.0 

5.0 

115.0 

9.0 

111.0 

14.0 

107.0 

18.9 

103,0 

24.0 

99.0 

29.7 

95.0 

TUNNEL  SPEED 


^hr-ft-°F^ 


36.14 
38.26 
38.45 
39 . 70 
41.00 
41.90 


36.8 


91.0 


41.81 


TABLE  19 


RESULTS  OF  INFRARED  SYSTEM  CALCULATIONS,  HIGH  TUNNEL  SPEED 


Run : 

#8 

Tunnel  Air  Velocity: 

306.0 

ft 

sec 

Average  Biot  Number: 

0.0022 

Average  Uncertainty 

in  h:  ±1.83 

. 

btu 

- 5 - 

hr- 

-ft  -°F 

Time  (Seconds) 

Temperature  ( 

°F) 

h 

.  btu 
^hr-ft 

_  o  ■ 

0.0 

119.0 

4,9 

115.0 

37. 

14 

10.3 

111.0 

37. 

46 

14.7 

107.0 

37. 

45 

20.6 

103.0 

37. 

20 

25.9 

99.0 

38. 

85 

31.5 

95,0 

40. 

45 

38.2 

91.0 

41. 

32 

TABLE  20 


RESULTS  OF  TVC  VANE  CALCULATIONS 

Run:  #1 

Tunnel  Air  Velocity:  214.3 

^  sec 

Average  Biot  Number: 

Average  Uncertainty  in  h: 

Time  (Seconds)  Temperature  ( ° F )  h 


0.0041 


±  3.23 


btu 


hr-ft  -°F 


0.0 

99.0 

5.5 

95.0 

38.52 

10.8 

91.0 

41;41 

17.0 

87.0 

41.89 

23.5 

83.0 

43.19 

31.6 

79.0 

43.31 

42.2 

75.0 

42 . 50 

53 


TABLE  21 

RESULTS  OF  TVC  VANE  CALCULATIONS 


Run : 


Tunnel. Air  Velocity: 

Average  Biot  Number: 

Average  Uncertainty  in  h : 


#2 


220.5 


ft 

sec 


0.0053 


±  4.41 


btu 


hr-ft  -“F 


Time  (Seconds)  Temperature  (°F) 


h  ( 


btu 


hr-ft-°F 


0.0 

97.0 

4.0 

93.0 

62.33 

9.0 

89.0 

59.10 

16.4 

85.0 

52.33 

23.2 

81.0 

53.63 

33.0 

77.0 

52.01 

47.3 

73.0 

49.17 

54 


TABLE  22 


RESULTS  OF  TVC  VANE  CALCULATIONS 


Run: 


#3 


Tunnel  Air  Velocity: 


220.5 


ft 

sec 


Average  Biot  Number : 


0.0042 


Average  Uncertainty  in  h:  ±3.90 


btu 


hr-ft  -°F 


Time  (Seconds ) 


Temperature  (°F) 


,  btu 
'hr-ft-°F 


LOW  SPEED 


HIGH  SPEED 


Figure  6.  Graph  of  Wedge  h  vs  Time,  High  Speed  Flow 


Flow  Direction  Relative  to  Picture 

t  t  I  t  I  t 

C 


E  TEMPERATURE  VS  TIME 
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